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Abstract 

We investigate a competition of tendencies towards ferromagnetic and incommensurate order in 
two-dimensional fermionic systems within functional renormalization group technique using tem- 
perature as a scale parameter. We assume that the Fermi surface (FS) is substantially curved, lies 
in the vicinity of van Hove singularity points and perform an account of the self-energy corrections. 
It is shown that the character of magnetic fluctuations is strongly asymmetric with respect to the 
Fermi level position relatively to van Hove singularity (VHS). For the Fermi level above VHS we 
find at low temperatures dominant incommensurate magnetic fluctuations, while below VHS level 
we find indications for ferromagnetic ground state. In agreement with the Mermin-Wagner the- 
orem, at finite temperatures and in small magnetic fields we obtain rather small magnetization, 
which appears to be a power-law function of magnetic field. It is found that the FS curvature is 
slightly increased by correlation effects, and the renormalized bandwidth decreases at sufficiently 
low temperatures. 

PACS numbers: 
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I. INTRODUCTION 



During last two decades, the problem of magnetic fluctuations in itinerant-electron com- 
pounds attracts substantial interest in connection with physics of layered systems. In 
copper-oxide high-temperature superconductors the importance of magnetic fluctuations 
is seen from the proximity of antiferromagnetically ordered region to the region of supercon- 
ducting state with pronounced incommensurate magnetic fluctuationsi»2.. Another impres- 
sive example of peculiar magnetic properties of low-dimensional systems is provided by the 
family of layered strontium ruthenates, which are pure metallic and incidentally have a very 
complex magnetic behavior. 

Neutron scattering- of the one-layer paramagnetic metal Sr2Ru04 (which at low tem- 
peratures T < 1.5 K becomes an unconventional, most likely p-wave, superconductor-i^) 
reveals an incommensurate character of magnetic fluctuation spectrum, with the largest 
contribution corresponding to the wave vector ~ (2/3)(7r,7r) and smaller one with the 
wave vector Q; ~ (l/3)(7r, tt). It was argued^ that incommensurate short-wave magnetic 
fluctuations characterized by the wavevector are expected to originate from a and fi 
bands (produced by dr^^ an dy^ orbitals and generating quasi-one-dimensional sheets of FS), 
while the long-wave magnetic fluctuations with the wavevector originate most likely 
from the 7 band, produced by d^y orbitals and generating quasi-two-dimensional sheet of 
FS. Moreover, magnetic properties of this compound can be controlled by the substitution 
effect: being doped by lantanium, Sr2_yLaj^Ru04 exhibits enhancement of magnetic suscep- 
tibility. This enhancement is likely related to the position of the Fermi level of 7 band: 
The fit of the tight-binding spectrum to experimentally determined Fermi surfaces allows 
to expect the onset of the ferromagnetic order at y = 0.27, related to the elevation of the 
Fermi level of the 7 band towards a van Hove singularity (VHS). However, no ferromagnetic 
order was found for this doping (Ref. 13), despite the VHS at the Fermi level, revealed by 
ARPES^. 

Bilayer ruthenate compound Sr3Ru207 also provides a number of interesting physical 
properties, in particular incommensurate magnetic fluctuations^ and metamagnetism^^; the 
latter was assumed to be related to the proximity of Fermi level to VHS^iii^. The effect of 



substitution of Sr by Ca in (Sri„a;Ca2:)3Ru207 was studied in Ref. 



13l . Alhtough the Wilson 



ratio approaches ~ 700 for x = 0.2, and the system becomes nearly ferromagnetic in the 
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temperature interval (3 < T < 10 K), further lowering of temperature does not result in 
the long-range ferromagnetic ordering, but instead forces the system to freeze into a spin- 
glass state. Therefore, an important problem related to the magnetic properties of metallic 
two-dimensional compounds is revealing the conditions, for which the ferromagnetic order, 
which is observed to be suppressed experimentally, is stable. 

From the theoretical point of view, magnetic properties of layered systems are closely 
related to their electronic structure. It is justified by ARPES experiments and first princi- 
ple calculations of above-mentioned compounds that the electronic spectrum is quasi-two- 
dimensional and can be fit using nearest (t) and next-nearest neighbour {t') hoppings. To 
investigate an interplay of different magnetic states in itinerant layered systems it is conve- 
nient to consider the two-dimensional (2D) one-band Hubbard model. Despite the fact that 
this model is strongly simplified as compared to the situation in real compounds, it catches 
an important dependence of ground-state ordering on the form of the electronic spectrum, 
in particular on the t'/t ratio and band filling. On the other hand, the electronic spectrum 
containing nearest and next-nearest neighbor hoppings retains VHS which is often present 
in the density of states (DOS) of real compounds. 

A large enough value of DOS at the Fermi level, occuring due to VHS, leads to ferro- 
magnetic ground state according to the Stoner theory^^. This approach, however, neglects 
magnetic fluctuations, in particular, of incommensurate type. A generalization of the Stoner 
theory to treat the spin spiral instability and Neel antiferromagnetism reveals that in a ma- 
jor part of the phase diagram the spiral phase is even more preferable than the ferromagnetic 
one^^. In particular, the phases with the wave vectors Q = {Q, Q) and Q = (Q, tt) were 
found to compete strongly with the ferromagnetic phasei^. The critical value of the Hubbard 
interaction Uc for the ferromagnetic instability is smallest in the region with the Fermi level 
below VHS, where the competition of ferromagnetic and (Q, vr) phase is observed. On the 
other hand, in the vicinity and well above VHS the ferromagnetic phase is to a great extent 
suppressed by the incommensurate (Q, Q) phase. 

To investigate the stability of ferromagnetic order with respect to long-wave magnetic 
fluctuations in itinerant 2D electronic systems at moderate Coulomb interaction, the spin- 
fermion model was recently used^"^. The static character of magnetic fluctuations in 2D 
systems with magnetically ordered ground state allowed to determine the boundaries of fer- 
romagnetic region and the most important magnetic instability, competing with ferromag- 



3 



netism. It was argued that the region of stabihty of ferromagnetism shrinks substantially, 
especially at VH filling, due to competition with incommensurate magnetic fluctuations, in 
qualitative agreement with mean-field^^ and variational results^. 

In the limit of large Coulomb repulsion the stability of saturated ferromagnetic ground 
state with respect to other magnetic structures was investigated starting from 1960s. In 
particular, Nagaoka has found that saturated ferromagnetic state is unstable with respect 
to the spin-wave excitation with the wave vector Q = (vr, vr), which implies that it competes 
with commensurate Neel antiferromagnetic order This result was later supplemented with 
variational principle results indicating that at least in part of the phase diagram the wave 
vector of the corresponding magnetic instability is incommensurate^^. 

Dynamical mean field theory (DMFT) provides a tool for unbiased account of the local 
electronic correlation effects and determination of the favorable type of ground-state mag- 
netic order for arbitrary interaction strength. Within DMFT the magnetic phase diagram 
for the Hubbard model was investigated on a Beth^^ and square*^ lattices with nearest 
and next-nearest neighbor hoppings. For small t' away from half filling and sufficiently 
large Coulomb interaction ferromagnetic, paramagnetic and incommensurate magnetic (the 
region where commensurate solutions of the DMFT equations do not converge) phases were 
found. At large enough t' < t/2 the ferromagnetism is found to be dominating in a broad 
region of parameters^i. 

The non-local electronic correlations for saturated ferromagnetic state can be treat in 
particular in the Kanamori's T-matrix approximation (which is exact in the limit of small 



carrier concentration), applied to this problem in Ref. |23|. However, this approach does 
not account for contributions from the other channels of electron interaction, which are 
important at not too small concentrations. The possibility of incommensurate magnetic 
ordering is also missed by this approximation. 

To study the interplay of magnetic and electronic properties and their effect on the pos- 
sibility of ferromagnetic instability, we focus here on the functional remormalization group 
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for 



(fRG) technique, which is a powerful tool for treatment of correlation effects (see Ref. 
the introduction in the application to the Hubbard model). In the pioneer study of Ref. 
the use of the Polchinski's form of fRG technique allowed to construct the phase diagram 
of the nearest-neighbor hopping Hubbard model {fi — T plane). The momentum cut-off 
scheme of the Wick-ordered version of fRG equations was employed to obtain the ground- 
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state instability type at small t' /t in Ref. |26l Later the temperature^"— and Hubbard 
interaction^'' were used as flow parameters in one-particle irreducible (IPI) investigations of 
the t — t' Hubbard model. Main result of these studies is that the small nearest-neighbor- 
hopping {f /t < 0.2) favors the antiferromagnetic instability at VH filling, moderate one 
(0.2 < t'/t < 0.35), the rf-wave supeconducting instability, rather large t'/t > 0.35 corre- 
sponds to ferromagnetic instability. Away from van Hove band filling the competition of 
antiferromagnetic and superconducting instabilities is obtained at small t'/ 1, while at larger 
values t'/t > 0.35 the competition of ferromagnetic instability and p-wave superconductivity 
is observed2I'2^'22. 

The essential shortcoming of the above reviewed approaches'""—*^ is that these do not 
consider the self-energy corrections to the electronic Green function. Hence, important non- 
trivial effects of renormalization of the electronic spectrum in the vicinity of magnetic phase 
transition (or in the regime with strong magnetic fluctuations) can be missed, provided that 
the Fermi level is near VHS. On the other hand, spin-dependent self-energy corrections 
provide the mechanism of the response to magnetic field, which is crucial in the context of 
invesigation of ferromagnetic instability. Note that only antiferromagnetic and supercon- 
ducting phases were considered previously within the combination of fRG and mean-field 
approach (Ref. l3ll ) and within the fRG approach in the symmetry broken phased. Consider- 
ing ferromagnetic instability poses a problem of scale-dependent FS, which was theoretically 
elaborated only in the self-adjusting Polchinski and Wick-ordered schemes, proposed in Ref. 
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In this paper we present a study of the evolution of magnetic and electronic properties 
with decreasing temperature within IPI fRG in zero and small finite magnetic field. We 
treat accurately the FS problem, including movement of projecting points and momentum 
dependence of the self-energy. The plan of the paper is the following. In Sec. |TT] we 
introduce the model and review earlier approaches to ferromagnetic instability at weak and 
intermediate coupling. In Sec. HH] we introduce the details of our novel fRG approach to 
take into account the fluctuation effects, retaining the electronic self-energy. In Sec. HVl we 
present and discuss the numerical results. In Sec. |V]we present the conclusions. 
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II. THE MODEL AND EARLIER APPROACHES 



We consider the Hubbard model with the action 

where /3 = 1/T is inverse temperature, is the number of lattice sites, /i is a magnetic field 
directed along z axis, measured in units of Bohr magneton, ek is an electronic spectrum, /i is 
the chemical potential, U is the Hubbard on-site interaction, 6 is the Kronecker 5-symbol. 
The sums in Eq. ([T]) are taken over 4-vectors k = (iz/„, k), where Un = 7i{2n + 1)T are the 
fermionic Matsubara frequencies (n G Z). We consider electronic dispersion on the square 
lattice 

ek = — 2t(cos + cos k^) + At' {cos k^ cos ky + 1), (2) 

where t,t' > are hopping parameters. Such a form of the spectrum corresponds to VHS 
at zero energy. We fix the ratio t'/t = 0.45. 

To gain insight into possible types of magnetic order, we consider the magnetic suscepti- 
bility Xq(w) in paramagnetic state^^, 

Xq'M=0q'M-t^, (3) 

where (pqioo) is the susceptibility, irreducible in particle-hole channel. The paramagnetic 
state is unstable with respect to the formation of incommensurate state with the wave 
vector Q, provided that the maximum of irreducible susceptibility 0q(O) is at q = Q and 
Xq(0) diverges. A critical interaction for stability of such an incommensurate magnetic state 
Uc = 0Q^(O). For the sake of simplicity let us consider the random phase approximation 
(RPA) where 0q(w) coincides with the non-interacting spin susceptibility 

Xq(i^n) = E G'^(i^n)Gtq(i^^n + iu^n) , (4) 

where G^liun) = {iVn — Ck + /^)~^ is the non-interacting electronic Green function. In 
Fig. [T] we present zero-temperature momentum profile of static non-interacting magnetic 
susceptibility Xql^^ = 0)- Below VH filling the competition between incommensurate mag- 
netic structure with large wave vector and a variety of incommensurate instabilities with 
small magnetic wave vectors is observed (see Fig. [1^). On the other hand, above VH fill- 
ing, the ferromagnetic instability competes with the long-wave incommensurate magnetic 

6 



Xoiq), t'=0.45t, iio=-0.05t 



Xo(q), t'=0.45t, fio=0.05t 












FIG. 1: Momentum dependence of the static magnetic susceptibility of free electrons x^i^ = 0); 
t' = 0A5t: a) Fermi level below VH filling (/x = -0.05t); b) Fermi level above VH filling = 0.05t) 



instability™(see Fig. [TJd). This consideration shows that for large enough t' the tendency to 
incommensurate magnetic ordering originates from the geometry of the FS. 

The competition of ferromagnetism and incommensurate magnetic ordering beyonc 
mean-field approximation was discussed within the quasistatic approach (QSA) in Ref. [17. 



Under the assumption of the ground state to be ferromagnetically ordered, it was found 
that magnetic fluctuations can reconstruct the momentum dependence of the susceptibility 
only at finite electronic interaction, stabilizing the ferromagnetic ground state with respect 
to incommensurate magnetic fluctuations. Above the critical value of the interaction, the 
maximum of 0q(O) is reached at q = 0. The corresponding result for the boundary sepa- 
rating the ferromagnetic and {Q, Q) incommensurate instabilities is presented below in Fig. 
El 

The effort of direct account of the electron correlation effects was performed in 
temperature-fiow fRG study21i2^ = 0-fRG). While the problem of ferromagnetism forma- 
tion was investigated in zero magnetic field, the self-energy renormalization was neglected, 
which makes the FS scale independent and the chemical potential not renormalized. Formal 
consequences of this assumption are discussed in Sec. IIIIDi 
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III. FUNCTIONAL RENORMALIZATION GROUP WITH SELF ENERGY 
CORRECTIONS 



Present study is an extension of the study of Ref. l28|. Contrary to previous approaches, 
we (i) do not neglect the self-energy corrections S to the electronic Green's function G and 
account for their momentum dependence, which results in the renormalization of spectrum 
parameters t and t' and moving FS, (ii) partially take into account momentum dependence 
of the vertex inside the patches, which allows to account for the scale dependence of the 
Fermi surfaces corresponding to different spin projections a = ±1. 



A. fRG equations 



We use the IPI fRG equations^^ in the truncation of Ref. 
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r = r4(G.G,).r 



PP 



r4(G.G,)*r 



ph 



r4(G.G,)*r 



phi 



(5) 



(6) 



where F is the IPI 4-vertex function (we refer to it as a vertex below), pp, ph, and phi 
denote particle-particle and two independent particle-hole channels, S is the self-energy, 
S = —G'j.dGQl/ds is the single-scale propagator, Go.s and Gg are appropriately rescaled 
non-interacting and interacting Green's functions, momentum arguments and spin indices 
are omitted for brevity. The asteriscs denote the summation over internal 4-momentum 
and spin indices corresponding to specified channel, dots denote the derivatives, taken with 
respect to the scale parameter s. 

If FM order parameter or magnetic field are directed along z axis, the SU(2) symmetry is 
broken. The remaining axial symmetry results in the self-energy which is diagonal in spin 
indices a, a' 

^aa'ik) = J^kaSaa', (7) 

and the vertex function components are nonzero provided that di + (T2 = o"3 + o"4- We 



consider the vertex component Ta-j^„^-aia2 
be obtained using exact relation 



(for brevity we denote it as Fg-^o-j), the others can 



r(Tio-2;cr3(T4 (^1) ^3? ^4) — To-i 



0"2 0"i;(T3(T4 



{k2, h; ks, = F 



(T2CTi;(T4(T3 



{k2, ki, ki, k3 



(8) 
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Another useful property, which is a consequence of the symmetry with respect to the com- 
bination of time-reversal transformation and a ^ —a transformation, is 



('■lO'2;cr3(T4 



(/ci, k2\ ki) — ro.30-4;(Ti(T2(^3, ^4; ki, k2)- 



(9) 



Neglecting the frequency dependence of the vertex and self-energy the fRG equations can 
be written in the following explicit form: 

rai(T2(ki,k2;k3,k4) = 

= {1-S^,^J2)^Y, Taia^lki, ka; p, ki+k2-p)£PP^^(p, ki+k2-p)r,,,,(p, ki+kg-p; kg, k4) 
p 

- ^ Xl^<^i<^'^^i'P'^3,ki -k3 + p)/:P^(p,ki -k3 + p)r^^2(ki - ks + p, ks; p, k4) 



+ ^<Tia2^ X] r^,<,(ki, k2 - ka + p; k4, ^)Cf^{p, kg - kg + p)!^^^!?, ^2] kg - ka + p, ka)] 

-(l-'^^aiaa)^ X] r<,,^2(ki, k2-ka+p; p, \^a)D'^^^^{p, k2-ka+p)r^i<,2(P, k2; ka, ka-ka+p) 
p 

(10) 



Ska = ^ r<^^'('^' P; ^' P) [/P^' + (2^P^' - SpcT')/pa'] - A XI r^^^'C^' P; ^' P)/p<T' - Ska/2, 
per' pir' 

(11) 

where 



^ka/L + ^k'a'/kV , (Ska " A)/ka + (^k'a' - A)/k'a' 



£^ka + £^k'a' 



+ 



£ka + ^k'a' 

(/ka + h'a' - l)(Ska + Sk'a' - 2/i) 
(£ka + £k'a')^ 



, (12) 



^^;'(k,k') = 



^kafLa ~ ^k'a'/k'a' , (^ka " A)/ka " (^k'a' " A)/k'a' (/ka " /k'a')(Ska " ^k'a') 



£ka ~ £k'a' 



+ 



£ka ~ £k'a' 



We have introduced the renormalized electronic spectrum 



£ka — e-k — iJ, — ha + Ek^., 



(Ska — £k'a')^ 

(13) 



(14) 



and fpa = /(£pa) = 1(1 - tanh(^£p^/2)). 
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Below we consider temperature cutoff^^ with s = log(t/T). The single-scale propagator 
has a form 

Sk. = -G,. y dT) ^^^^ 

where Gka = (iJ^n — £ka)~^- In the present study we choose temperature independent bare 
chemical potential {dfi/dT = 0), although the temperature dependence of /i can be adjusted, 
e. g., to keep the number of particles fixed^S. The equations ffTOl) and ffTTl) are basic for IPI 
fRG approach since they can, in principle, determine the renormalization of the electronic 
spectrum and interaction at any temperature. This system of equations should be sup- 
plemented by initial conditions at s = — oo(T = +oo). In the infinite-temperature limit 
correlation effects are absent, and one obtains for the self-energy Sko- = U/2 and interaction 
r^i^2(fci, k2; h, ki) = U{1 - ^aiaa)- 

The fRG equations (fTOj) . (ITT]) form the system of ordinary functional differential equations. 
To solve them numerically, one has to introduce some approximation procedure to reduce the 
considered system to a finite system of ordinary equations. Below we present the procedure 
of numerical solution in details. 

B. The self energy ansatz and Fermi surface 

The first step is the use of common patching scheme^^ to avoid dealing with system of 
functional equations. However, we do not neglect the momentum dependence of the self- 
energy (and vertices) inside the patches, in particular Eq. ffTTj) is fully employed. We assume 
that the self-energy has the form 

Sko- = — 2(5to-(cos + cos ky) + A5t'^ cos k^ cos ky + T,„. (16) 

where the parameters St„,6t'^ and are determined as follows. While solving numerically 
the system of fRG equations, at each step we calculate the values of Sko- on two sets of 
(o"-dependent) projecting points (PPs, see for details Appendix). For considered k G PPs 
the linear regression procedure is used, 

Sko- — )■ — 25to-(cos kx + cos ky) + Adt'^ cos kx cos fcy + Sg-. (17) 

In this way we determine the flow of unknown quantities 5ta, 6t'^ and So- in Eqs. (ITOl) . (ITTl) . 
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Such a choice efficiently reduces the number of variables and retains VHS points of the 
renormalized spectrum. Note that does not depend on k and renormalizes the chemical 
potential /i; 5tcr and 5t'^ contribute to the change of momentum dependence of the spectrum 
{5tfj corresponds to the bandwidth renormalization). 

Therefore, it is convenient to represent the renormalized spectrum (IT^ in the form 

£^k<7 = -2teff,<7(cOS K + COS K) + 4t[,g ,^(C0S cos ky + l) - fieS,a, (18) 

where 

Atefr,a = fl-T.a + 4:6t'^ + ha, teS,a = t + 6t^ t'^g^^ = t' + 6t'^. (19) 

The scale-dependent effective chemical potential fies,a results in Fermi surface, calculated 
with the renormalized spectrum parameters teff.o-, t'cs a- ^PPlyi^S present method one should 
be careful in determining the geometry of the FS. If t^g^/tefr.o- < 1/2, the bottom of the 
band is at the energy = — 4teff,(T + 8t'^s „ singly connected; however, if 

^eff o- Aeff.cr > 1/2 the bottom of the band is zero. In this case if fies,a > Wa-, the FS is 
singly connected and the patching scheme of Ref. bsl can be used, while for /ies.o- < w^, the 
Fermi surface consists of two disconnected parts and the patching scheme should be chosen 
differently. 

Despite that the self-energy renormalization are included in the present study, we neglect 
incoherent contributions to the Green functions. Due to this renormalization, the actual 
Fermi level fies is determined by the combination of the bare spectrum and the self-energy 
parameters f|T9|) at the end of the flow. 



C. The vertex ansatz 



The vertex function F is represented by its values at the current FS. Since we always 
trace the renormalization of Ska- we have to take into account moving of the Fermi surface 
during fRG flow. PPs of the current FS are changing during the flow and the discrete 
(projected) vertex function derivative acquires an additional contribution corresponding to 
this movement, 

dV _ dV dV dkpp 
ds ds dkpp ds 

We denote symbolically the derivatives with respect to PPs (/cpp) as d/dkpp. 
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To take into account this momentum dependence of T we assume that, apart from the 
position of external legs in certain patches, the vertex function depends on momenta kj 
through the renormalized energies e^-^ and linearize the latter dependence. Let the momenta 
of external legs k'[,k2,k3 be on the current FS and consider vertex with all momenta kj 
belonging to the main set of PPs and three vertices with two momenta kj belonging to the 
main set of PPs (see Appendix), and one belonging to the auxiliary set. Therefore we have 
4 possibilities for the choice of ki,k2,k3, and obtain a system of four linear equations (we 
use ek^a, = 0) 

r..,(k5, ki; ki)) + J2 5^r..,(k5, ^; k^))^k,.. = r..,(ki, k2; kg), (21) 

i 

with unknown quantities T„„i(k1,'k2;k.l)) and (9jro-o-'(k'j;, k2; kg). Solving the system (|2T|) 
one determines the vertex on the current FS. 

D. Relation to previous approaches 

In zero magnetic field, neglecting the momentum dependence of the self-energy (i. e. with 
6t = 6t' = 0), one can demand the chemical potential to absorb all the corrections to the 
self-energy in Eq. ( fTOl) . which is possible since the electronic spectrum enters the Eq. ( JTOjl 
through ^ko- only. On the other hand, Eq. ffTTj) . being reformulated with /i' = yU — S, reduces 
to /i' = 0. Hence, in this case, /i is replaced by a constant The important consequence 
of this approximation is that FS is fixed within this ansatz, since it is determined solely by 
the scale-independent parameter fi'. 

If the field is non-zero, the self-energy corrections cannot be absorbed into the chemical 
potential even when the renormalization of the hopping parameters is neglected. This, in 
turn, results in moving of spin FSs during the flow, so that projected vertex acquires a 
contribution from this moving. As described in detail in Sect. IIII C\ we account for this 
contribution in the self-consistent numerical scheme of treatment of Eqs. (fTOl) . ( ITTj) . and 

(EH). 
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IV. TEMPERATURE DEPENDENCE OF THE RENORMALIZED PARAME- 
TERS 



In this Section we present and discuss numerical results of the present fRG approach 
accounting for the self-energy corrections in zero (Sec. IIV Ap and finite magnetic field (Sec. 
IIVBI) . Afterwards, in Sec. IIV CI we present our results for the phase diagram. 

A. Results in zero magnetic field 

In this subsection we consider the results of fRG calculations in the spin symmetric phase 
for U/t = 3 and 4. We choose the bare chemical potential /i in such a way that at the end 
of the flow the renormalized position of the Fermi level /les, determined by Eq. (1T9|) . lies in 
the vicinity of VHS. It is clear physically and verified numerically in our investigation that 
the magnetic response is suppressed for the Fermi level well separated from VHS. 

Starting from an infinite-temperature limit the position of the Fermi level first tends 
to increase from its high-temperature Hartree value fi — U/2 to almost low-temperature 
Hartree value Ji = fi — U j^^i^^ii p{e)de, p being the bare DOS. At even lower temperatures 
fi decreases due to correlation effects. The obtained low-temperature scale dependence of 
for different bare chemical potentials fi is shown in Fig. [2^ {U/t = 3) and Fig. {U/t = 4). 
We stop the flow, when the effective interaction becomes too large (> 2.5W, where W = 8t 
is bare bandwidth). The corresponding scale Smin yields the minimal temperature which is 
available within the flow, Tmin = ^ exp(— Smin)- In the following we parametrize the initial 
chemical potential of flow by the Hartree Fermi level fl. 

The cases U = 3t and U = 4t are somewhat different due to the absence of the region 
of saturation of fies{s) dependence for U = 3t, although the dependence fief[{s) becomes 
weak at the end of the flow. Below we consider the renormalized Fermi level fi^s in the 
saturation region (if it exists) as a representative parameter which characterizes the flow, 
since the renormalization of other parameters of the electronic spectrum is not too strong. 
The charge response is slightly suppressed in the vicinity of van Hove filling, for the case 
fieff < more than for fi^s > 0. Non-monotonous behavior of fics{s) slightly above VHS 
{fi = l.Ot) for U = it, which consequences are discussed below, is worthy of notice. 

In our scheme the effective chemical potential renormalization has substantial contribu- 
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Effective chemical potential, U=3t Effective chemical potential, U=4t 




s=ln(t/T) s=ln(t/T) 

FIG. 2: Temperature dependence of the renormalized Fermi level fi^s at /i = 0: a) U = 3t, b) 
U = At. Chemical potential values p. (see text) are shown by numbers near the plots. 

the spectrum parameter tg^/t, U=4t, h=0 the spectrum parameter tg^'/tg^, U=4t, h=0 




s=ln(t/T) s=ln(t/T) 

FIG. 3: Effective (renormalized) hopping parameters: a) tes/t, b) tgg/tcff for U = At. Numbers 
correspond to the value of chemical potential fi 

tion from t and t' renormalizations (see Eq. ( fT9l) ). In Fig. [3] we present for instance tes/t and 
t'es/tcs plots for U = it. While the bandwidth (tes/t) is somewhat reduced at the end of the 
flow well below VHS, indicating prominent correlation effects; for fi^s above VHS tcs/t first 
increases with decreasing temperature (which implies that correlation effects in this regime 
are not substantial) and then decreases at lower temperatures, which is related to enhance- 
ment of the vertex in the vicinity of magnetic instability. On the other hand, the ratio t'/t 
is not strongly renormalized (this means that the self-energy effects do not change substan- 
tially the curvative of the Fermi surface) and monotonously increases towards the value 1/2 
as /i increases, which corresponds to flattening of the electronic dispersion. The complete 
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U=3t 



U=4t 




1.7 1.8 1.9 

s=ln(t/T) 



FIG. 4: The scale profiles of T^*^^ (solid lines) and F^^^'^ (dashed lines), see text, in the case h = 
a) U = 3t, b) [/ = 4t for different choices of the chemical potential p, (shown by numbers) yielding 
Fermi level below VHS {ficS < 0) 



treatment of the case > lAt, where tes/^eff exceeds 1/2 during the flow, corresponding to 
a change of FS geometry (see for details Sec. IIIIBI) . requires special consideration and is 
beyond the scope of the present study. 

The type of leading magnetic instability can be inferred from the behavior of vertices. 
Let us consider the scale profiles of maximal vertex. We consider maximal total vertex T"^^^ 
(the maximum is taken over all possible combinations of momenta) and maximal exchange 
vertex F™**^'^ (the maximum is taken over all combinations of momenta with k2 = ks). The 
ferromagnetic instability is accompanied by coincidence of T^^^ and F^^^'^ in the vicinity of 
transition point where both the values diverge, reflecting an instability of the paramagnetic 
state with respect to zero-momentum collective spin excitations. Therefore, the criterion 
for commensurate magnetic fluctuations has the form 



AF, 



-\max.E 







(22) 



and yields an information about the type of leading instability. 

The scale dependences of T^'^'' and T'^^'^'^ for different /i are presented in Fig. S] {^^s < 0) 
and Fig. \5\ (yUeff > 0). Let us consider first the case of Fermi level below VHS, /icfr < 
0. For both the cases U = 3t (Fig. H^) and U = At (Fig. |Dd) F^'^'^ is diverging, and 
AFmax/r^^'^ vanishes or is very small at the end of the fiow. Therefore, magnetic fiuctuations 
are predominantly ferromagnetic at the end of the fiow. On the other hand, AFmax increases 
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FIG. 5: The same as in Fig. |4]for Fermi level being above VHS (/Xcff > 0) 

as /ieff approaches VHS: for both U = 3t and U = At AFmax is the largest in the case Ji = 0.8t. 
This means that in the vicinity of VHS the ferromagnetic fluctuations hardly dominate over 
incommensurate ones. 

For fics > (see Fig. [5]), the difference AF^ax is non-zero up to the lowest temperatures 
where Fmax diverges. Moreover, the ratio AFmax/r^i'*'^ increases as /i increases. We interpret 
this as that the incommensurate fluctuations are dominating over the ferromagnetic ones. In 
the case ft = l.Ot, where the temperature dependence of /ieff(s) is strongly non-monotonous 
(see Fig. Wp), which is possibly related with invalidity of ansatz f fT6|) . fl2T|) in this regime and 
causes the shift of the point of vertex diverging to substantially lower temperatures. This 
case should be considered more correcly in further elaborated studies. 

In the next subsection we introduce a small magnetic field to investigate magnetic prop- 
erties of the system. 



B. Magnetization in finite magnetic field 

In this subsection we supplement the picture in zero magnetic field considered above by 
the results for magnetic response in finite magnetic field. In case h > we have /Xefrt f^eSi 
due to spin dependence of the spectrum parameters tefr.o-, ^es o- ^o- and, to a small extent, 
due to the presence of magnetic field. At low scales (high temperatures) the spectrum 
parameters do not depend on the spin projection substantially and fi^s^ — fi^si ~ ^h, but 
with increasing the scale (lowering temperature) the strong spin dependence can be realized 
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FIG. 6: The logarithmic plots of the dependence m{h) and its fit to m cx h". The chemical 
potentials p, and fitted a are shown near the plots. 

which is a manifestation of exchange enhancement. 

The strengh of magnetic response is characterized by the magnetization m which can be 
easily calculated using current parametes of the electronic spectrum St^j, 6t'^, So- (we remind 
that the frequency dependence of the self-energy is neglected), 



Fig. [7] shows the comparison of the scale profiles of magnetization m for different bare 
chemical potentials. For /Xes above VHS we find slightly negative m at the end of the 
flow for both the cases U = 3t and U = 4t. This is possibly related to an influence of 
incommensurate magnetic fluctuations above VHS, as conjectured from the results in zero 
magnetic fleld. Below VHS, where dominating ferromagnetic fluctuations were observed in 
the absence of magnetic fleld, the magnetic response to magnetic fleld becomes positive and 
considerable, especially in the vicinity of VHS: the maximal value of m at the end of the 
flow increases, but the temperature of sharp increase of the magnetization becomes lower as 
/i increases. Note, that the absolute values of magnetization are rather small (m <C 2n) at 
the end of the flow. To verify the fulflllment of Mernin- Wagner theorem, we flt the data for 
magnetization taken at lowest temperatures of the flow to m oc /i" dependence (see Fig. [6]), 




(23) 



k 
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FIG. 7: The magnetization scale profiles in the field h = 10 for different chemical potential fi 
specified by numbers: a) U = 3t, h) U = At 

which gives a G (0.46, 0.83) depending on the chemical potential /i. In the vicinity of VHS 
level (/2 = l.Ot, l.lt) for U = 4t the magnetization is somewhat larger at the end of the flow, 
in particular a = 0.17 at /i = l.lt. Additionally, in this case magnetization tends to saturate 
at lowest temperatures with increasing the magnetic field (not shown). For these chemical 
potentials the renormalized Fermi level lies very near VHS and the Mermin- Wagner theorem 
may not be fulfilled to a good accuracy. On the other hand, these cases are on the borderline 
between the regions of strong ferromagnetic and incommensurate fluctuations. 

The temperature dependence of the exponent a is of particular interest, since it allows to 
investigate the magnetic properties of the system while entering the region of strong magnetic 
fluctuations. It is obvious that at high temperatures m oc h {a = 1), but at low temperatures 
this relation can be violated due to fluctuations. To determine the temperature of the 
crossover into the regime with strong ferromagnetic fluctuations, we adopt the criterion a = 



18 




s=ln(t/T) 



s=ln(t/T) 



2.05 




600 
500 
400 
f 300 
200 
100 




h=5 lO^^t 
h=1 10 ^ 
h=2 10,t 
h=4 lO'l 




0.3 
0.25 

0.2 
0.15 

0.1 
0.05 




I . 
h=5 lO^t 
h=1 10,t 
h=2 10,t 
h=4 lO'l 



H=0.4t 

b) 



1.7 1.75 1.8 1.85 1.9 1.95 

s=ln{t/T) 



s=ln(t/T) 



2 2.05 




2.05 2.1 




1 2.2 



ln(t/T) 



s=ln(t/T) 



FIG. 8: The scale profiles of ratios m/h (a,c,e) and m/h}^^ (b,d,f) for U = 4:t at different small 
magnetic fields: h/t = (5,10,20,40) • 10"^ and different bare chemical potentials /i, which values 
are shown explicitly 

1/3, dictated by the mean-field value of the critical exponent at the conventional magnetic 
phase transition. We denote the crossover scale as s* (the corresponding temperature is 
T* = texp(-s*)). 

Fig. [S] shows the comparative (with respect to the magnetic field) representation of m/h 
and m/h^^^ scale profiles for different positions of Fermi level below VHS in the case U = At. 
In cases /i = 0.4t, 0.6t the crossover scale is reached within the fiow and lies near its end, while 
in the case /2 = 0.9t it is beyond the fiow and is obtained by an extrapolation. Therefore, 
increasing /2 tends to decrease T* . We do not find a well-resolved crossover scale at t/ = 3t, 
but we note that the system is close to it at the end of the fiow in the case /i = 0.4t, where 
the magnetic fiuctuations tend to be commensurate. This means that at f/ = 3t crossover 
to ferromagnetic ordering can be realized only at extremely low temperatures which are not 
available within the present approach. For the case /Xeg > and both U = 3t and U = it 
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FIG. 9: Summary magnetic phase diagram in the vicinity of VH fihing in n^s — U plane, t' = OABt, 
obtained within the mean-field^(short-dashed hnes, = M ~ Un/2, where n is the electronic 
density), S = O-fRCi^ (solid lines), QSA— (long-dashed lines) and present SE fRG (symbols) 
approaches, see text. The boundary lines corresponding to S = 0-fRG and QSA approaches 
separate ferromagnetic and paramagnetic regions below van Hove filling and ferromagnetic and 
incommensurate magnetic phases above van Hove filling. 

we find a relatively weak response which becomes negative close to the end of the flow, as 
discussed above. 

The distribution of temperatures T^in and T* for different chemical potentials yields a 
natural connection between zero-field (Tmin) and finite-field (T*) results. We find that 
both Tinin and T* decrease as the system approaches VHS, which is possibly an effect of 
incommensurate magnetic fluctuations and temperature smearing of VHS (above van Hove 
filling DOS falls down faster than below it); below VHS we obtain T* < Tmin. However, 
in the case /i = l.lt we find T* > T^in, which is possibly connected to the abovediscussed 
violation of the Mermin-Wagner theorem in a close vicinity above VHS. 
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C. Phase diagram 



In this Section we summarize the results obtained in zero (Sect. IIV AO and finite f llVBp 
magnetic field and compare our results with the results of previous approaches. The phase 
diagram constructed in terms of renormalized Fermi level /igg (see discussion in Sec. IIV Ap 
and Coulomb interaction U is shown in Fig. |9l 

In the case fi^s < the large region on the phase diagram is found where the zero-field 
magnetic fiuctuations are predominantly commensurate. In finite magnetic field, the system 
exhibits some indications of ferromagnetic order in the ground state: at f/ = 3t a crossover 
from paramagnetic to ferromagnetic order is observed well beyond the fiow, while aX U = 4t 
the crossover to ferromagnetic state is well-resolved in the close vicinity of the end of the 
fiow. In the case fies > (excluding the region of near vicinity of Fermi level to VHS) we do 
not find commensurate magnetic fiuctuations to be dominating in zero magnetic field: they 
are competely replaced by incommensurate fiucuations with corresponding maximal vertex 
F^^^ being diverging. This conclusion agrees with the result in finite magnetic field, where 
the magnetization m becomes slightly negative at low temperatures. Therefore we conclude 
that ferromagnetic order is suppressed by an incommensurate magnetic fiuctuations if the 
Fermi level is above VHS. The case of fies very near but above VHS is worthy of special 
attention: in the case U = 3twe find incommensurate magnetic fiuctuations in zero magnetic 
field and no indication of ferromagnetic ordering in finite field. However, very near VHS 
at f/ = 4t we find almost commensurate magnetic fiuctuations in zero magnetic field and 
ferromagnetic-like behavior in finite magnetic field. At the same time, in this case fieS 
depends non-monotonously on fi and the magnetization on the value of magnetic field. 
These states are shown by "?" symbol in Fig. [91 

The mean-field and quasistatic approximation (MFA and QSA) boundary lines calculated 
in Refs. jl6j and il7\ are shown on the phase diagram in Fig. [9] for comparison (the effective 
chemical potential in MFA = Un/2 and in QSA it is determined from the condition 
for electronic density n(/iefr) = noifi), where no(yu) is the number of particles in the non- 
interacting model). The critical interaction for stability of ferromagnetism within the present 
approach Uc > St is somewhat larger then the QSA and the mean-field resuts {UdficS = 
0) ~ 2t), which naturally suggests that the account for the electronic correlations results in 
an enhancement of Uc- 
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The results of the present study partly agree with those calculated within fRG approach 
without self-energy effecta^^ (E = 0-fRG approach, see Fig. [9]), obtained by studying the 
temperature dependences of magnetic and superconducting susceptibilities in zero magnetic 
field. The lower threshold for ferromagnetism Uc obtained in the present study, appears to 
be finite contrary to S = 0-fRG. At the same time, away from VHS these two approaches 
qualitatively agree. 

The comparison of different approaches MFA, QSA, S = 0-fRG, SE fRG shows a step- 
by-step restriction of the size of ferromagnetic region. Ferromagnetism is practically absent 
for /iefT > within SE fRG, but not restricted with respect to previous approaches in the 
case yUeff < 0. Above VHS ferromagnetism is destroyed by well-resolved incommensurate 
magnetic fluctuations, while below VHS quantum commensurate fluctuations dominate. 

V. CONCLUSIONS 

In this study we present a fRG treatment of magnetic order in the Hubbard model, 
controlled by the Fermi level being in the vicinity of van Hove singularity (VHS) and on-site 
Coulomb interaction. The introduced version of fRG accounts for self-energy corrections, 
which implies a proper account of the Fermi surface moving as the temperature decreases. 
The following aspects of this problem are investigated: the renormalization of vertices in zero 
magnetic field, spin splitting in finite field, finite-temperature behavior and the electronic 
spectrum renormalization. 

We find that magnetic properties of the system are substantially asymmetric with re- 
spect to the Fermi level position relatively to VHS. In particular, ferromagnetic ordering 
is strongly suppressed above VHS due to competition with incommesurate magnetic fluc- 
tuations. Below VHS we find precurors of ferromagnetic ordering at low temperatures and 
not too small interaction (of order of half bandwidth). The temperature of the crossover to 
strong ferromagnetic fluctuations decreases as the Fermi level elevates towards VHS. 

Correlation effects do not change substantially the form of electronic spectrum which is 
characterized by t'/t. However, the renormalized ratio t'/t monotonously increases as the 
Fermi level rises towards VHS. When the Fermi level is at VHS we find that the ansatz 
used for the self-energy momentum dependence is too crude to catch delicate features of the 
renormalization. 
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The results of the present approach improve the results of the mean-field approximation^ 
and quasistatic approach™, and agree with previous fRG study^ for the position of the Fermi 
level away from VHS. The Mermin- Wagner theorem is also shown to be fulfilled to a good 
accuracy in the present approach in the most part of the phase diagram. The magnetic 
field dependence of magnetization demonstrates power-law behavior with the exponents 
a G (0.62,0.83). 

The observed region of ferromagnetism for the Fermi level above VHS is rather narrow. 
This explains qualitatively the magnetic behavior of one- and two-layer ruthenates (which 
have, after corresponding particle-hole transformation, /igfr > 0). The non-Fermi-liquid 
behavior of doped Sr2Ru04'' can be possibly related to the non-monotonous temperature 
dependence of the Fermi level near VHS in the present approach, which neglects non- 
quasiparticle contributions. These contributions may therefore be important for the Fermi 
level near VHS (see, e. g. Ref. A, which is the subject of future study. The results obtained 
demonstrate an important role of many-electron renormalizations of electron spectrum, in 
particular of chemical potential, in the presence of VHS. This fact can be crucial for the 
criterium and properties of weak itinerant ferromagnetism. Another important point is 
that peculiarities of ferromagnetic ordering related to the presence of VHS cannot be even 
qualitatively captured by any elaborate method based on DOS consideration only. 

The obtained finite-temperature picture, in particular, the relation of the character of 
magnetic fluctuations to the position of the Fermi level and absence of the long-range order 
(which is stated by the Mermin-Warger theorem) should be possibly supplemented by a 
zero-temperature study of this problem. Investigation of the effect of finite (but not too 
small) magnetic field and possibility of the metamagnetic transitions is also of interest. 
Another unsolved problem is a unified description of weak and strong ferromagnets, which 
can receive new insights from weak-coupling investigations of the Hubbard model. 

We are grateful to W. Metzner and H. Yamase for stimulating discussions. This work was 
supported in part by Partnership Program of the Max-Planck Society, the project 'Quantum 
Physics of Condensed Matter' of Presidium of Russian Academy of Science and President 
Program of Scientific Schools Support 4711.2010.2. 
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Appendix A: Projecting points 



To parametrize the momentum dependence of the self-energy and vertex, we introduce 
a set of points in the Brillouin zone (below we refer to them as projecting points, PPs), so 
that the values of functions at these points represent the function. In prevous studies^I this 
set of the points was chosen as an intersection of lines of constant angles (located at the 
center of patches) with FS which was assumed to be fixed. Since we account for momentum 
dependence of the self-energy and vertex more accurately, we supplement this set of PPs 
at FS by a corresponding set in the vicinity of FS by choosing additional points belonging 
to "shifted" FS and call it auxiliary set of PPs. However, since moving of FS is continuous 
process, we apply discrete Runge-Kutt procedure (RKP) to solve the system of differential 
equations numerically, which causes some complications in definition of PP. 

At any discrete step of RKP we introduce the PPs of the following types: (i) main PPs, 
which belong to FS, determined by the chemical potential fi, and the self-energy Sko- at the 
beginning of the step, (ii) auxiliary PPs which are detemined analogously to main PPs, with 
fi being shifted by 6fi„ for different spin projections, and is determined by typical shift 
of momentum-independent part of self-energy So- at previous step, (iii) current PPs which 
are determined by current Sko- within the step. Note, that at the beginning of the step main 
PPs coincide with current ones. Introduction of the auxiliary PPs set is needed since, due 
to definition of main PPs set, the functions chosen to represent momentum dependence of E 
(cos kx + cos ky, cos k^ cos ky, 1) are linearly dependent at the constant energy set (see Section 
nil Bp . which does not allow to use them for linear regression of S-derivatives. Therefore, 
an additional set of PPs is needed; the calculation of the vertex near current FS allows 
to expand the vertex linearly beyond simple projecting ansatz (see detail explaination in 
Section nil C|) . which was used in earlier studies. 

RKP makes its step using the vertex and self-energy derivatives calculated at the inter- 
mediate (current) value of argument s at fixed (during the step) sets of main and auxiliary 
PPs. From this we extract derivatives 6tcr, 6t'^ and So-, see Section IIII B[ However, the 
right-hand side of fRG equations contains current self-energy function and vertex function, 
projected on the current FS. RKP allows to calculate this current self-energy function and 
use it to determine current FS and corresponding current PPs set. Then we apply the pro- 
cedure of Sect. IIII CI accounting for the infiuence of the FS moving on the current value of 
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projected vertex. 
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